Charge-transfer excitations in molecular donor-acceptor complexes within the 

many-body Bethe-Salpeter approach 
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We study within the perturbative many-body GW and Bethe-Salpeter approach the low lying sin- 
glet charge-transfer excitations in molecular donor-acceptor complexes associating benzene, naphta- 
lene and anthracene derivatives with the tetracyanoethylene acceptor. Our calculations demonstrate 
that such techniques can reproduce the experimental data with a mean average error of 0.1-0.15 eV 
for the present set of dimers, in excellent agreement with the best time-dependent density functional 
studies with optimized range-separated functionals. The present results pave the way to the study 
of photoinduced charge transfer processes in photovoltaic devices with a parameter-free ab initio 
approach showing equivalent accuracy for finite and extended systems. 



PACS numbers: 31.15.A, 33.20.Kf, 71.35.-y 

A large set of studies aiming at clarifying the processes 
assisting the separation of the photogcncrated electron- 
hole pairs in organic solar cells revealed the importance of 
charge transfer excitations at donor /acceptor interfaces i 
From a theoretical point of view, the fundamental issue of 
correctly modelling, at the quantum mechanical level, the 
electron-hole interaction, is a long standing difficulty. In 
particular, while standard density functional theory^ can 
accurately describe excitations with a strong spatial over- 
lap between the involved occupied and unoccupied states, 
the interaction of weakly overlapping electron-hole pairs 
requires the development of novel functionals with a fine 
tuning of the local and non-local exchange and correla- 
tion contributions^ a strategy followed e.g. in the recent 
developments around the range-separated functionals^ 

In this work, we study within the many-body GW 
and Bethe-Salpeter (BSE) perturbation framework the 
charge-transfer excitation energies of ten donor/acceptor 
molecular complexes, associating tetracyanoethylene 
(TCNE) with benzene, naphtalene, anthracene and their 
derivatives. We show that these techniques provide re- 
sults which are in excellent agreement with available ex- 
perimental data. The accuracy is equivalent to that of 
the best time-dependent density functional theory based 
on optimized range-separated functionals. We emphasize 
in particular the importance of self-consistency when per- 
forming GW-BSE calculations starting from DFT calcu- 
lations with (semi)local functionals. 

The choice of such donor-acceptor complexes com- 
bining TCNE and acene derivatives is dictated by the 
availability of gas phase experiments in the case of ben- 
zene, toluene, o-xylene and naphtalene donors^ comple- 
mented by experiments in solution for anthracene and its 
derivatives^ Such experiments provide invaluable refer- 
ence data for the study of the merits of the various the- 
oretical approaches, even though experiments performed 
in solution complicate the comparison with theoretical 
calculations performed on the isolated complexes. 

We adopt the geometries obtained by Stein and 
coworkers at the DFT-B3LYP level in their TDDFT 
study of the same systems using an optimized 





FIG. 1: (Color online) Isocontour representation of the (a) 
HOMO and (b) LUMO states of the anthracene-TCNE com- 
plex. The contours have been set to 20% of the maximum 
value of the wavefunction. Blue and red contours indicate dif- 
ferent signs of the wavefunction. The grey, white and brown 
atoms are carbon, hydrogen and nitrogen, respectively. The 
ir/n* character of the wavefunctions is apparent. 



range-separated functional labeled BNL(7=0.3) in the 
following^ To illustrate the typical geometry of such sys- 
tems, and highlight their donor-acceptor character, we 
plot in Fig. 1 an isocontour representation of the high- 
est occupied (HOMO) and lowest unoccupied (LUMO) 
molecular orbital of the anthracene-TCNE complex, 
showing the clear localization of the HOMO (LUMO) 
on the anthracene donor (TCNE acceptor). 

As compared to the TDDFT approach, the GW-BSE 
formalism proceeds in two steps. First, the quasiparti- 
cle HOMO-LUMO gap, namely the difference between 
the ionization energy and electronic affinity, is calculated 
within the GW approximation^ - — Traditionally, such 
calculations are performed non-self-consistcntly, namely 
as a "single-shot" correction to the starting DFT-LDA 
Kohn-Sham eigenstates, an approach labeled Go Wo in 
the following. While this approach has been shown since 
the mid-eighties to yield excellent results for most ex- 
tended semiconducting systems,— several studies^— ev- 
idenced that in the case of small molecular systems, the 
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GoWo(LDA) approximation yields too small ionization 
potential and HOMO-LUMO gaps. It was shown that 
such a problem could be cured either by starting from 
Hartree-Fock eigenstatesji2r— or by performing a sim- 
ple self-consistent update of the quasiparticle energies 
when building the self-energy operator , 15 ' 16 an approach 
labeled GW in what follows. 

In a second step, the neutral (optical) excitation en- 
ergies can be obtained as the eigenvalues of the Bethc- 
Salpetcr H e ~ h Hamiltoniar>i2r— which, in the |</>f<^ > 
product basis of the non-interacting unoccupied \<pl > 
and occupied > single-particle states, is composed of 
three terms as follows: 



(a) Gas phase experiment 



(b) Experiments in solvant 
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with (ef P ,(t)^ h ) the best single-particle eigenstatcs in 
the absence of electron- hole interaction. W(r,r') and 
V (r, r') are the (statically) screened and bare Coulomb 
potential respectively Following the above discussion, 
the ef P can be taken to be the G W (LDA) or GW 
quasiparticle energies. Clearly, a too small HOMO- 
LUMO gap in H ag will result in too small neutral ex- 
citation energies. 

Our GqWq-BSE and GW-BSE calculations are per- 
formed with the Fiesta codei^ which exploits an even- 
tempered auxiliary basis representation of the two-point 
operators such as V c and W, while the input DFT- 
LDA cigenstates are taken to be those of the Siesta 
DFT code22 using a large triple-zeta plus double polariza- 
tion (TZ2P) basis. In the case of the smaller complexes 
with benzene, toluene and o-xylene, these calculations 
are double-checked using the planewave basis set2S. 

We compile in Fig. 2 and Table I our calculated excita- 
tion energies within the Go Wo (LD A) -BSE (yellow down 
triangle) and GIT-BSE (black diamonds) approaches, 
that we compare to the experimental data (red cir- 
cles) and to the TDDFT results of Stein and coworkers^ 
with the B3LYP (blue up triangle) and range-separated 
BNL(7=0.3) (green squares) functionals. Following the 
analysis by Stein and coworkers^ experimental data ob- 
tained in solution (light red circles) have been shifted 
upwards by a 0.32 cV constant that approximately ac- 
counts for the solvent induced bathochromic shift. As 
shown in the Table, our results as obtained with the Fi- 
esta code and with the planewave Yambo approach (in 
parenthesis), are in good agreement within 0.1 eV accu- 
racy, showing that the present results are independent of 
the technicalities of the various implementations. 

The mean absolute error (MAE) as compared to ex- 
periments is found within the GW-BSE approach to 
be 0.10 eV for "gas phase" benzene, toluene, o-xylcne 
and naphtalene complexes, and 0.13 eV averaging over 
all systems, including the corrected data from the ex- 
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FIG. 2: Experimental and theoretical first singlet excitation 
energies in eV. The signification of the symbols is given in (a). 
Experimental data are from (a) Ref. [f for gas phase experi- 
ments on benzene, toluene, o-xylene and naphtalene donors, 
and (b) Ref. @ for experiments in solution on anthracene and 
its derivatives. Following Stein and coworkers in Ref. 0, a 
constant 0.32 eV shift has been added to the experimental 
data in (b) to (tentatively) mimic the bathochromic shift. 
The B3LYP and BNL data are the theoretical TDDFT re- 
sults from Ref. 0. Figures (a) and (b) are not on the same 
energy scale. 



donor 


G W -BSE 


GW-BSE 


BNL 


Exp. 


benzene 


3.03 (3.10) 


3.58 (3.63) 


3.8 


3.59 


toluene 


2.52 (2.67) 


3.27 (3.37) 


3.4 


3.36 


o-xylene 


2.23 (2.28) 


2.89 (3.00) 


3.0 


3.15 


naphtalene 


1.96 


2.55 


2.7 


2.60 


MAE 


0.74 (0.68) 


0.10 (0.07) 


0.12 




substituent 


anthracene derivatives 


none 


1.66 


2.17 


2.14 


(2.05) 


9-carbo-methoxy 


1.64 


2.05 


2.16 


(2.16) 


9-cyano 


1.89 


2.32 


2.35 


(2.33) 


9-methyl 


1.21 


1.99 


2.03 


(1.87) 


9,10-dimethyl 


2.00 


2.21 


2.09 


(1.76) 


9-formyl 


1.89 


2.32 


2.27 


(2.22) 


MAE 


0.43 


0.15 


0.11 





TABLE I: Experimental and theoretical optical gap (eV) for 
donor-TCNE complexes. G W -BSE and GW-BSE results 
in parenthesis have been obtained with a planewave basis im- 
plementation (see text). MAE is the mean absolute error. 
The BNL column reports the TDDFT results with optimized 
range-separated BNL(7=0.3) functionals from Ref. [j. Exper- 
imental results are from Ref. he The experimental results in 
parenthesis were obtained in solution and a 0.32 eV constant 
energy has been added (see text). 
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pcrimcnts in solution. An interesting and satisfactory 
result is the good agreement with the previous TDDFT 
calculations based on range-separated functional, with 
a localization parameter 7=0.3 (see Ref. 0) optimized 
to reproduce the correct quasiparticle band gap. In the 
two cases where theory and experiment start showing a 
severe disagreement, namely with methyl and dimethyl 
anthracene derivatives (the two lowest experimental ex- 
citation energy in Fig. 2b), both TDDFT-BNL and GW- 
BSE results agree in overestimating the extrapolated ex- 
perimental data. This may indicate that the estimation 
of the solvation effect is incorrectly accounted for by the 
rigid 0.32 eV shift. Finally, as emphasized above, an 
important observation is that the GqWq-BSE approach 
provides significantly too small excitation energy, as a re- 
sult of too small a GoWo-LDA quasiparticle band gap<2! 
As expected, the GW-BSE results are much better than 
those obtained with the TDDFT-B3LYP approach^ 

In conclusion, we have shown that GW-BSE first- 
principles many-body perturbation theory provides 



charge transfer excitation energies in excellent agreement 
with experiment, with a mean average error of the or- 
der of 0.1-0.2 eV. The present results are of equivalent 
accuracy as compared to the best range-separated func- 
tionals TDDFT calculations, with an approach showing 
equivalent merits for finite size and extended systems 
without any adjustable parameter. In fact the screened 
Coulomb potential W(r, r') through which the holes and 
electrons are interacting is a non-local operator with a 
decay in finite size or extended systems automatically ad- 
justed through the evaluation of the dielectric properties 
of the system. The analysis of such interaction in donor- 
acceptor complexes may help in the future to better un- 
derstand the features of the best exchange-correlation 
kernels within the TDDFT framework. 
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